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Abstract 

The modulational instability (MI) in the class of NLS equations is 
discussed using a statistical approach (SAMI). A kinetic equation for 
the two-point correlation function is studied in a linear approximation, 
and an integral stability equation is found. The modulational insta- 
bility is associated with a positive imaginary part of the frequency. 
The integral equation is solved for different types of initial distribu- 
tions (5- function, Lorentzian) and the results are compared with those 
obtained using a deterministic approach (DAMI). The differences be- 
tween MI of the normal NLS equation and derivative NLS equations 
is emphasized. 
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1. Introduction 

The modulational instability (MI - also known as Benjamin- Feir in- 
stability) is a well known phenomenon encountered, in quite general 
conditions, almost each time a quasi- monochromatic wave propagates 
in a weakly nonlinear medium (Dodd et al 1982, Remoissenet 1999, 
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Scott 2003, Benjamin and Feir 1967, Bespalov et al 1966 and Abdul- 
laev et al 2002). It is responsible for the generation of robust, localized 
excitations, the solitons in completely integrable systems, and solitary 
waves in non-integrable ones. 

Two different ways to discuss this phenomenon are possible. The 
first one is a deterministic approach (DAMI) in which the amplitude 
of a plane wave, with an amplitude dependent dispersion relation (a 
Stokes wave), is slowly modulated 

t) = a(l + eA(x, t)) exp (i(kx - ut)) (1) 

This approach is well known and can be found in any textbook on non- 
linear physics (see Dodd et al 1982, Remoissenet 1999). The second 
approach is a statistical one (S AMI) , and its aim is to see the influence 
of the statistical properties of the medium on the MI phenomenon. It 
is a less used approach, but in the last ten years was intensively ap- 
plied to investigate different physical situations (Alber 1978, Janssen 
1983, Fedele et al. 2000, Hall et al 2002, Lisak et al 2002, Fedele et 
al 2002, Onorato et al 2003), ranging from hydrodynamics (Janssen 
1983, Fedele et al. 2000, Hall et al 2002, Lisak et al 2002, Fedele et al 
2002), to plasma physics (Fedele 2002), large amplitude wave-envelope 
propagation in nonlinear media, high-energy accelerators, condensed 
matter. Recently it was extended also to discrete systems (Visinescu 
et al 2003, Grecu D. et al 2005) and to systems of coupled nonlinear 
equations (Manakov's system) (Grecu D. et al 2004 and 2005). 

In this approach a kinetic equation for the two-point correlation 
function 

p(xi,x 2 ) = <*(xi)**(x 2 )) (2) 

is written down and a linear stability analysis near an initial distribu- 
tion 

Po(xi,x 2 ) = po(\x 1 -x 2 \) (3) 

is performed. The initial state is assumed homogenous and conse- 
quently po depends only on \x± — x 2 \. In the previous expression of 
p(xi,x 2 ) by (. . .) we understand a statistical average over the distri- 
bution function of the medium. 

The aim of this paper is to investigate MI, especially using the 
statistical approach of the following equation 

.99 ld 2 9 XIT|2t . lTl2 0tf . d|^| 2 T n 
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It is an extended NLS equation, and for different choises of the real pa- 
rameters A, fi, v several completely integrable equations are recovered. 
If fi = v = 0, we get the well known NLS equation 

,_ + __ +W4 , = (5) 

while for A = 0, fi = v one obtains the completely integrable derivative 
NLS-I equation (dNLS I) (Kaup 1978, Nakamura et al. 1980, Mio et 
al 1976, Mj0lhus 1976) 

.m IdH . d / 1ti2t \ 

and for A = v = we get another integrable derivative NLS II equation 
(dNLS II) (Chen, Lee 1979) 

.dV ld 2 ^ . lTl2 d# n 

Before starting the SAMI study of eq. ((3J), let us briefly recall the 
results of a DAMI analysis. Introducing in Q and keeping only 
the linear terms in e a linear evolution equation for A(x, t) is found. 
Looking for plane wave solutions of it 

A(x,t) = ai e i{Qx - m) + a* 2 e-' l[ - Qx ' n * t] 

the compatibility condition of the homogenous system leads us to 



n-Q[k + (n + i/)|a| 2 ] = iQy (A - u\a\ 2 - jjik)\a\ 2 - ^ (8) 

an the instability is related to ImSl > 0. For A / and denoting 
Ai = A — Hal 2 we obtain 



M = Qy(Ai - l^k)\a\ 2 - ^ (9) 

and one sees that the instability regions depend on the sign of fi. In 
a (Q 2 , k) representation these are given in Fig. IT"! (arbitrary units) 

It is interesting to consider the instability regions for the two inte- 
grable dNLS equations. For dNLS I an instability arises only if fi < 0. 
We have 

lmn = Qj\iu\\a\ 2 (k-\fi\\a\ 2 )-^- (10) 
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Figure 1: Instability regions (A / 0); (a) fj, < 0, (b) fj, > 0. 



and in a (Q 2 , k) representation the instability occurs only if k > 
|/i||a| 2 . For dNLS II the same condition, [i < 0, must be fulfilled 
and 



ImSl = Q]J \[i\\a\ 2 k — 



(11) 



The instability regions in (Q 2 , k) representation (arbitrary units) are 
presented in Fig. 2. 





(b) 



Figure 2: Instability regions of the integrable dNLS equations: (a) dNLS I, (b) 
dNLS II 

These results have to be compared with the similar result obtained 
for NLS equation ©, which writes 



ImQ = Q\ X\a\ 2 - ^ 



(12) 
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We see that in the NLS case the instability region is independent 
on the wave vector k of the carrying wave, and restricted to the long 
wave-length region (Q 2 < 4A|a| 2 ), while for all the situations discussed 
above they are ^-dependent, and more over they are not restricted 
to the long wavelength region (exception is the case (b) of Fig 1). 
Actually for the dNLS equations the long wavelength region is either 
forbidden, or highly unfavourable for the development of the MI. 

In the next section the kinetic equation for the two-point correla- 
tion function is obtained and its stability is studied in the linear 
approximation. An integral stability equation will be obtained and 
solved for different initial distributions in section 3 ((^-distribution, 
Lorentzian). Few conclusions and comments are given in the last sec- 
tion. Preliminary results have been presented in Grecu, A.T. 2005. 



2. Kinetic equation. Linear stability 
analysis 



The kinetic equation for the two-point correlation function (j2J, is ob- 
tained easily in the following way: 

1. write the equation @ for x = x\ and multiply by $?*(x 2 ); 

2. write the complex conjugated of (J1J) for x = X2 and multiply by 



3. add the two equations and take the statistical average; 

4. finally use a Gaussian approximation to decouple the four-point 
correlation functions. 

As an illustration we give below two examples of such decouplings 



{^{xi)^>*(xi)^(xi)^*(x 2 )) = 2a 2 (x 1 )p{x 1 ,x 2 ) 




= a- 



2 (xi)t, — p{xi,x 2 ) + p{xi,x 2 ) lim ^ — p(x 1 ,x 2 ) 




where 



a 



■ 2 (xi) = <tf(xi)tt*(xi)> 



(13) 
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The kinetic equation found in this way is 



.dp 1 ( cP_ d_ 

1 di + 2 \OxJ ~ dx 2 



p + 2X a 2 (xi) — a?(x2) 



+ v) 
+2ivp 



— d — 

a 2 (xi)- h a 2 {x 2 )^— 

ax\ 0x2 



d 8 
nm — h nm — — 

Xl-*X\ Qx\ xi^x 2 0x2 



9 v 9 

iim — h nm — — 

X 2 ^Xl ijXl X!^X 2 OX2 



P 
P 

p* = 



The next step is to introduce the new set of variables 
x = x\ — X2 relative coordinate 
X = - (xi + X2) center of mass coordinate 



(14) 



(15) 



and make a Fourier transform with respect to the relative coordinate. 
This procedure is known as the Wigner-Moyal transform (Wigner 
1932, Moyal 1949) (also known as Klimontovich's statistical aver- 
age method, Toda et al 1995), and is very useful to deal with non- 
homogenous evolution equations. The Fourier transform is defined 
as 

' 1 * dxe-' lkx p{x,X,t) (16) 



Using the definition ljlM|) of a 2 , it is easily seen that 



a 2 {X,t) 



dkF(k,X,t) 



(17) 



After straightforward calculations the Fourier transform of (|14j) is 
given by 



OF OF A ™ 
+ k— + 4A^ 



■1)> d 2 i +1 a 2 (X)d 2 i +1 F(k,X) 



Of, 



dX 



j=0 



(2j + l)!2 2 J+ 1 OX 2 ^ 1 



Ok 2 i +1 



0a 2 (X) 9 



a 2 (X) 



F(k,X) 



(18) 



OX ' V ~'" J ' ~ ox 
™ (-1)3 2 ^(X) 2 iF{k,X) 

+ KP + V)1^ 2j[2 2j QX 2j Q X Q k 2j 
3=0 J 

_ 2(m +„)jr / 2j+ Y (x)a2j+1[fc 5 fc ' x)] = 



^(2j + l)!2 2 i+ 1 OX 2 'J+ l 



Ok 2 ^ 
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In obtaining (JT5|) we have expanded a 2 {x\) = a 2 (X + 1) and a 2 (x2) = 
a 2 (X — |) in power series around the point X, and we have used the 
relation 

Also we used 

and similar expressions for the other limits which appear in l|14|l. 
For a linear stability analysis we consider 

F(k,X,t) = f(k) +eT(k,X,t) (19) 

and 

a? (X, i) = a§ + eof (X, t) (20) 

where 

r+oo r+oo 

a 2 = dkf(k), al(X,t)= dkF(k,X,t) (21) 

J — oo J — oo 

Here f{k) is the Fourier transform of pod^l), which, due to the ho- 
mogenous assumption of the initial state, depends only on k, and 
moreover is an even function of k. 

The linearized kinetic equation in the (k,X) representation is given 

by 

W + [k + (M + »)<® m + (M + ^ + 2,/ j 

_ 9 r + if (- 1 ) 3 ' a 2j+ W(^)g 2j+1 (fc/(fc)) _ 
^ ^ (2j + i) ! dx 2 ^ 1 

The last step is to look for a plane wave solution of (|22|) 

F(k,X,t) = 9l {k)S QX -^ + 5 *(fc)e- i ( QX - n *') (23) 
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where 



Gi 



/+oo 
gi (k)dk, » = 1,2 
-oo 



Introducing (|23|) into the following equation is found 

(k - u)g\ + (p + v)fG\ + 2vfG 2 
2A 



+ 





H + v 



Gi 



*(*+-)-*(*-- 



d = 



(24) 



(25) 



Here we denoted w = % — (fi + u)oq and h(k) = kf(k). Also we used 



- (2j + 1) ! V 2 / dk^ +l 



and a similar relation for /i(fc). 

In (|25jl we have to divide by (A; — u>) and integrate over k. Intro- 
ducing the notations 



I 



/(*) 



dk 



J 



K 



_oo k — UJ 

+OG /(* + %)- f(k 



— oo 



k — uj 



dk 



(26) 



+oo h(k + §) - h{k - 

-oo & — UJ 



dk 



the following homogenous algebric linear equation in G\ and G2 is 
obtained 

Gi + (/x + + 2i/JG 2 + ^ JGi - ^^KG X = (27) 

A second equation is found starting from the complex conjugated of 
namely 

G 2 + (// + u)IG 2 + 2vIGx + ^ JG 2 - f^^KG 2 = (28) 
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The compatibility condition for the system ()27 [) .(|28 j) leads us to the 
following integral stability equations 

l + ^-u)I + ^J-^K = (A) 

or 

l + ^ + 3u)I+^J-^K = (B) 

In the next section these will be solved for different initial conditions 
/(*)• 



3. Solution of the integral stability equa- 
tion 

The first distribution function we shall discuss is a (5-spectrum 

f(k) = ^S(k) (29) 

It doesn't describe a realistic situation as it coresponds to a constant 
two-point correlation function in the initial state, po = Oq = const.. 
But the calculations are easily done and can be considered as a limit 
case to which more realistic situations can be compared. The integrals 
/, J, K are immediatly calculated 

r _j$ tJ __&^ iK _ (30) 

Then both integral equations (A) and (B) become an algebraic equa- 
tion of third order in u) 

— — Q 2 —O 2 

uj 3 - pa 2 oo 2 + (2\al - 2-)u + paf^ = (31) 

where p = a*a = A* ~~ v f° r A-equation and p = \xb = A 4 + 3f for 
B-equation. With u> = y + |/* a o it is reduced to the canonical form 

y 3 +py + q = (32) 

where 

P = 2X^-\{p4) 2 -^ 
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2 -- 



_2 l/r-._2\2 , Q 



Xa 2 



(my + 



(33) 



Complex solutions of (jSU) are obtained if the discriminant A = (|) 3 + 
(|) 2 is positive. The sign of A is controlled by the sign of p. For A 
we have p > if Q 2 /4 < 2Aoq — |(/iao) 2 , corresponding to the long 
wavelength region. With the notations 



IP Q 
— = -tana, tano? 
27 2 



a 

tan — 
2 



we have 



Imw 



sin 2(/? 

But if p < we can still have A > 0. Writing 



91 
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2Aa^ 



(Aa2) 2 + 3Q 2 



the condition A > leads us to the inequality 



^ 2 - —(fla 2 ) 2 



Q 3 - \ia\Q 2, < fial 
In this case (A > and p < 0) the usual notations are 



p° q 

-— = — sin a, tanw 
27 2 ' ^ 



tan ^ 



and one obtaines 



Imw 



-p cot 2tp 



The instability remains in the long wavelength region. 
When A = the situation is completely changed. Then 



P 



Q 



l 



^ + §(/^o) 2 



2 S 



Q 2 

4 
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(34) 
(35) 

(36) 

(37) 

4331) 
(ESI) 



(38) 



and A > can be realized only if Q > |/2aQ. This means that the 
instability develops for higher wave vectors Q, a result in complete 
agreement to the conclusions of the DAMI analysis of the derivative 
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NLS equations. Using the same notations (34') we get the same result 
(35'), and the requirement A > reduces to 

i) 2 >-(!) 3 

an inequality from which the limit value of ^- can be calculated. 
A more realistic situation is a Lorentzian distribution 



7T k + <J 

which corresponds to an exponential decaying two-point correlation 
function in the initial state 

Pq{x) =^e~ ax x > (41) 

The integrals I, J, K are easily calculated in the fc-complex plane 
assuming u complex with Imw > 0. We obtain 



a, 



UJ 



J = (42) 



K _ -fropQ 



where Co = to + ia, and the instability condition Imw > becomes for 
the new variable Co, ImcD > a. 

To have a point of reference, let us recall the result for the usual 
NLS equation (// = v = 0) (Visinescu et al 2003, Grecu D. et al (2004, 
2005), Grecu A.T. 2005) 



Imw = W2Aag - ^ -a (43) 



We get \m.uj > if 



Q < 2^2Xal-a 2 (44) 

which shows that the MI is restricted to the long wave-length region 
and is strongly dependent on the correlation length in the initial state; 
if a 2 > 2Xal, the initial state is stable to small modulations. 



11 



When Q42JI are introduced in (A) and (B), we get again an algebraic 
equation of third order in to, but now with complex coefficients. It can 
be reduced to the canonical form with 

Q 2 



P 



X + \ va ) a l ~ 3W) - — 



(45) 



(^o) 2 + 



91 
4 



where 7i has the same meaning as in the case of the 5-spectrum, and 
V = [i+v. As mentioned before the MI region coresponds to Imu) > a. 

Instead of working on the general case let us consider only the 
situations of the dNLS equations. We shall start with the dNLS-II 
equation (A = v = 0). Then both equations (A) and (B) reduce to a 
single equation 



T2 1 



1 - /uog- + 



a 







LO 



(46) 



It is convenient to measure 00, aandQ in units of |/u.|a§. If e = sgn//, 
the following algebraic equation of third order is found: 



lO ' 



iea ] lo + e 



Q 2 







(47) 



It has a certain symmetry having solutions of the form uj = eui r + iu)j. 
As we are interested in the imaginary part of these solutions, the 
result is independent on the sign of fi. The equation (|47|) is solved 
numerically. Measuring Q in units of (//a^) 2 , one has ImQ = Q (1m lo — 
a). A three dimensional plot , ImQ function of a and Q 2 /4 is given 
in Fig. ED As one sees, the neighbourhood of Q = is excluded from 
the instability region, which extends on a finite domain of Q and a. 
One can find a limit value for a {ac — 0.194) and the MI exists only 
for a < ac- 

Now let us consider the case of dNLS-I equation (A = 0, fj, = v). 
The equation (A) becomes 



l + 2i-^L 
which is easily solved giving 



a) 2 







(48) 



Imd; = —= 
V2 



t) 2+4 « 2 "t 



1/2 



(49) 
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Figure 3: Instability region for dNLS I. Three dimensional plot of ImQ versus 
Q 2 /4 and a; Two different orientantions. 

and the MI region is determined by 

, a 2 < nal (50) 

This instability region is located near the origin in the plane (Q, a). 
Another instability region results from (B), which writes 

cD 3 - 4aD 2 - ^ - 2i£^ u + = (51) 

where as in the case of dNLS-II, we measure u), a and Q in units of 
|/i|aQ. This equation has the same form as (|46|) . The instability region 
is obtained numerically and is represented in Fig. 4. Again the Q ~ 
region is excluded and one finds a critical value a c ~ 0.441 for a, which 
determines the superior limit of the MI region. 

4. Conclusions 

Several conclusions result from this analysis. For the extended NLS 
equation @ the instability region in a DAMI analysis is dependent on 
the wave-vector k of the carrying wave. In the case of derivative NLS 
equations the long wave-length region is either forbidden, or highly 
unfavourable for the development of MI. 



Q 2 < 4 
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Figure 4: Instability region for dNLS II. Three dimensional plot of Imfi versus 
Q 2 /4 and a; Two different orientantions. 



In the SAMI instead of the ^-dependence the instability region is 
now dependent on the correlation length in the initial 2-point corre- 
lation function. Although we got explicit results only for a 5— and 
Lorentzian distribution function we conclude that the MI is possible 
only for long-range correlations in the initial state. In the case of a 
5-function initial distribution, in a limit situation corresponding to an 
infinite correlation length, the results are similar with those obtained 
from DAMI analysis; for the dNLS equations again the long wave- 
length region is highly unfavourable. The same conclusions can be 
drawn for the more realistic case of a Lorentzian initial distribution 
function, namely that the MI region is restricted to a finite region 
of (Q, a) plane containing the axis a = 0. A critical value a c is 
determined for both dNLS equations and again the Q ~ region is 
unfavourable for the instability development. 

Another realistic initial condition is a Gaussian distribution function 



Calculations done only for NLS eq. (see Alber 1978) reveal the same 
conclusions, namely the MI is possible only for long-range correla- 
tions in the initial state (a small) and in the long wave-length region. 
The result is more complicated, being expressed through some special 
functions (complex integral function u(z), Abramowitz and Stegun 



Such calculations for dNLS equations are not yet performed. 



f(k) = a 2 



1 




e 2, 



1965). 
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The influence of the statisticai properties of the medium on the 
development of large amplitude waves, is a problem of special interest 
in hydrodynamics, especially for the generation of freak waves in the 
ocean. In a linear theory there is no interaction between ocean waves. 
A focusing phenomenon of wave energy may occur only when con- 
structing interference takes place. The situation is competely differ- 
ent when nonlinear wave- wave interaction is present. Then a wave can 
borrow energy from its neighbours and due to this extra-focusing phe- 
nomenon, waves with very high amplitude can be generated. In this 
process the knowledge of the initial conditions of the oceanic waves 
is essential. Accordingly, many experimental works have been con- 
ducted. For a long time the Joint North Sea Wave Project (JON- 
SWAP) has studied the power spectrum of the oceanic waves (Komen 
et al 1994), resulting a non-Gaussian expression. It was used to calcu- 
late the probability to generate large amplitude oceanic waves (ana- 
lytically and numerically) in given realistic conditions (Janssen 1983, 
Onorato et al 2001, 2003, 2004). 

Calculations using other initial distribution functions, reflecting 
different physical situations, would be of real interest for sure, even in 
the case of derivative NLS equations. 
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